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ABSTRACT  cracks  emanating  from  the  fin  tip  along  its  own  axis  of 


Computational  analysis  and  two-dimensional  tensile 
tests  on  single  motor  grain  fins  suggest  that  cracks  in 
fin  tips  are  most  likely  to  originate  at  the  coalescence  of 
a  fin  end  tip  radius,  with  a  small  radius  from  the  side  of 
the  fin.  Some  manufacturers  have  also  noticed  defects 
formed  during  casting  at  the  fin  tip  on  the  fin  axis. 
The  following  is  an  experimental  investigation  utilizing 
frozen  stress  photoelastic  models  of  an  existing  motor 
grain  geometry  in  order  to  clarify  stress  intensity 
factor  (SIF)  values  and  crack  growth  paths  for  cracks 
emanating  from  the  two  above-  noted  potential  critical 
loci.  Comparisons  between  results  from  cracks  grown 
from  the  two  loci  will  be  made,  suggesting  interesting 
conclusions. 

INTRODUCTION 

Two-dimensional  thermal  shrinkage  tests  on  pho¬ 
toelastic  models  and  two-dimensional  tensile  tests  on 
single  fin  sections  of  motor  grain  have  suggested  that, 
for  a  specific  fin  geometry  consisting  of  a  small  edge  fin 
radius  coalescing  into  a  large  fin  tip  radius,  the  criti¬ 
cal  locus  for  crack  initiation  lies  at  the  point  of  con¬ 
fluence  for  the  two  radii  at  the  fin  tip.  On  the  other 
hand,  some  motor  grain  manufacturers  have  reported 


symmetry  as  a  result  of  voids  collecting  there  during 
the  casting  process.  One  aspect  of  the  problem  which 
has  not  been  explored  is  how  cracks  would  grow  from 
these  two  critical  points  in  a  three-dimensional  model 
and  the  corresponding  stress  intensity  factors. 

The  first  author  and  his  colleagues  have  achieved 
some  success  (Smith,  1990),  (Smith  and  Wang,  1992), 
(Smith,  Wang  and  Liu,  1995),  (Smith,  Liu  and  Wang, 
1995)  combining  the  frozen  stress  photoelastic  method 
(Appendix  A)  together  with  a  two-parameter  method 
(Smith  and  Kobayashi,  1993)  for  analyzing  local  stress 
fields  around  crack  tips  for  (SIF)  determination  in 
generic  and  specific  rocket  motor  geometries.  The 
present  paper  describes  an  attempt  to  use  this  method 
to  study  the  crack  growth  and  SIF  values  for  cracks 
emanating  from  the  two  critical  loci  described  above 
in  three-dimensional  models  under  internal  pressure. 

TEST  MODELS  AND  PROCEDURES 

All  tests  models  were  cast  at 
Measurements  Group  (Raleigh,  NC)  in  molds  provided 
by  Virginia  Tech  using  PLM-4BR  stress  freezing  ma¬ 
terial  manufactured  by  Measurements  Group.  The 
model  dimensions  are  pictured  in  Fig.  1.  The  starter 


cracks  were  made  by  .first  drilling  a  small  hole  opposite 
the  fin  in  which  the  crack  was  to  be  located,  sliding  a 
shaft  with  a  tip  blade  into  the  hole,  positioning  the 
blade  at  the  critical  point  on  the  fin  surface  and  then 
striking  the  shaft  with  a  hammer.  The  starter  cracks 
then  emanated  from  the  blade  tip  into  the  material  as 
“natural”  or  real  cracks.  For  the  cracks  located  at  the 
point  of  confluence  of  the  small  fin  edge  radius  with 
the  larger  fin  tip  surface  radius,  a  shorter  tip  blade  was 
used  opposite  the  main  blade  to  balance  the  impact 
load  on  the  shaft  (Fig.  2).  This  shorter  blade  caused 
some  slight  damage  and  sometimes  a  small  crack  which 
was  patched  with  PMC-1  cement.  Of  course,  cracks 
emanating  from  the  critical  locus  on  the  axis  of  sym¬ 
metry  required  no  balance  blade  and  were  more  easily 
produced.  After  the  starter  cracks  were  made,  the 
drilled  holes  were  plugged  and  sealed. 

It  is  important  that  crack  growth  be  achieved 
from  these  starter  cracks  under  pressure  in  order  to 
establish  the  “pressure”  crack  growth  path  as  opposed 
to  the  starter  crack  generated  by  wedge  forces. 

The  test  setup  is  pictured  in  Fig.  3.  The  pressure 
was  monitored  only  outside  the  model,  as  indicated, 
and  maintained  its  value  when  end  caps  leaked  or  when 
a  crack  grew  through  the  model  wall.  In  addition  to 
the  internal  pressure,  an  end  load  was  applied  to  the 
model  to  approximately  simulate  the  effect  of  an  open 
end  as  the  motor  grain  may  experience  in  service. 

Four  identical  models  have  been  tested  to  date. 
Model  #1  was  utilized  solely  to  determine  the  stress 
fringe  order  at  the  two  critical  loci  without  cracks  (Fig. 
4).  Since  the  test  material  showed  a  slight  deteriora¬ 
tion  in  the  surface  layer,  graphs  of  fringe  order  versus 
path  along  “squeezed”  fringes  were  extrapolated  to  the 


boundary  and  showed  little  difference  in  fringe  value 
at  the  two  critical  points  (Fig.  5). 

Models  #2,  #3  and  #4  each  contained  two 
cracks;  in  each  case  in  fins  adjacent  to  an  uncracked 
fin  and  one  with  a  plugged  hole  in  order  to  allow  inde¬ 
pendent  analysis  at  each  crack  (Fig.  6).  After  stress 
freezing  and  removal  of  thin  slices  mutually  orthogonal 
to  the  local  crack  surface  and  its  border,  the  slices  were 
analyzed  by  a  two-parameter  model,  briefly  described 
in  Appendix  B. 

RESULTS 

As  noted  above,  two  types  of  cracks  were  used  in 
the  experiments.  Cracks  which  were  symmetric  with 
respect  to  both  load  and  geometry  were  initiated  at 
the  fin  tip  along  the  axis  of  the  fin.  These  cracks 
were  planar,  semi-elliptic  cracks  which  maintained 
these  characteristics  as  they  grew.  SIF  values  were 
determined  at  maximum  depth,  along  or  near  to  the  fin 
surface,  and  exhibited  pure  Mode  I  SIF  values.  Two  of 
these  cracks  penetrated  the  outer  surface  of  a  model 
and  readily  grew  at  smaller  load  levels  than  the  off- 
axis  cracks  which  emanated  at  the  locus  of  confluence 
of  the  two  fin  tip  radii.  These  off-axis  cracks  were 
generally  non-planar  and  were  extremely  difficult  to 
grow  even  at  elevated  pressures.  This  was  because  the 
tips  of  these  cracks  were  under  mixed  mode  load  and 
had  to  turn  a  significant  amount  in  order  to  grow.  In 
fact,  with  one  exception,  none  of  these  cracks  reached 
a  purely  Mode  I  state  during  testing. 

The  test  data  and  results  are  presented  in  a  nor¬ 
malized  form  in  Table  I.  Since  stress  freezing  pressure 
was  maintained  after  breakthrough,,  data  did  not  ap¬ 
pear  to  be  adversely  affected.  However,  fringe  observa¬ 
tions  on  Model  1  and  SIF  data  on  one  end  of  the  long 


Table  I 


Loads1 


P  -  3 1 1 .4  N  Long  off*axis  crack 

Pmax  =  0.1 45  MPa  a/c  =  0.66  a/t  =  0.26 


p^Tia*  i  =  1,  2 


Crack  Description2 

depth  (d) 

surface  (s) 

top 

bottom 

Model  2 

Long  off*axis  crack 

CO 

O 

(N 

II 

cC 

2.42 

2.90 

pSf=  0.07  MPa 


P  =  311.4N 


Long  symmetric  crack 
a/c  =  0.9  a/t-  1.0 

Model  3 

Short  symmetric  crack 


F2  =  0 
penetrated 


fin  tip  surface 
I  (top) 

/  |S  Tnc 


Tt . 

fin  axis 


Pmax  ~  0.103  MPa  a/c  =  0.5  a/t  =  0.41 
Psf=  0.033  MPa 

Long  symmetric  crack 
a/c  =  0.68  a/t  =1.0 


penetrated 


P  =  88.97  N 


Model  4 

Right  off-axis 


Pmax  =  0.049  MPa  a/c  =  0.81  a/t  =  0.25 
psf=  0.035  MPa 

Left  off-axis 
a/c  =  0.78  a/t  =  0.23 


Fj  =  1.51+ 
F2  =  0.47+ 

Fj  =  1.90 
F2  =  0.48 


(bottom) 

ns  =  near 
surface 

t  =  distance 
from  fm  tip 
to  outer 
surface  along 
fin  axis, 
(37.08  mm) 


+  Opposite  crack  surfaces  were  not  uniformly  separated  through  the  slice  and  this  led  to  local 
fringe  distortions. 

Notations: 

1  -  P  —  axial  compressive  load 

Pmax  =  maximum  internal  pressure  to  grow  crack 
psf-  stress  freezing  pressure 

2  -  a  =  crack  length 

c  =  half  length  of  crack  in  fin  tip  surface 

3  -  V5  =  approximation  of  elliptic  integral  of  second  kind 

(  V'65 

£=1  +  1.464  -  -<1 

W  c 

(  V  65 

£=1  +  1.464  -  ->1 

All  flaws  were  characterized  as  semi-elliptic  flaws  of  depth  a  and  length  2c.  However,  off-axis 
cracks  were  neither  perfectly  semi-elliptic  nor  planar. 


symmetric  crack  in  Model  2,  together  with  surface  val¬ 
ues  for  the  right  off-axis  crack  in  Model  4,  suggested  a 
small  skin  effect  in  the  material,  so  these  values  were 
omitted.  In  order  to  escape  this  effect,  values  of  Fi 
were  obtained  20°  off  of  the  surface  or  near  the  fin 
surface  “ns”  in  Model  4  and  appeared  quite  reason¬ 
able.  All  of  the  cracks  grew  more  along  the  axis  of  the 
model  than  through  the  thickness. 

The  long  off-axis  crack  in  Model  2  was  the  only 
such  crack  to  grow  sufficiently  to  escape  the  shear 
mode.  It  is  pictured  in  Fig.  7. 

SUMMARY 

Based  on  the  above  studies  to  date,  it  has  been 
found  that: 

i.  Cracks  emanating  from  a  fin  tip  on  the  axis 
of  the  fin  exhibit  pure  Mode  I  and  grow  more 
readily  than  those  emanating  from  the  point  of 
confluence  of  tip  radii. 

ii.  Cracks  emanating  from  the  confluence  of  tip  radii 
exhibit  mixed  mode  conditions  around  the  crack 
front  and  maintain  a  mixed  mode  condition  as  the 
crack  front  turns,  exhibiting  strong  resistance  to 
further  crack  growth. 

iii.  Normalized  stress  intensities  based  upon  a  planar 
semi-elliptic  shape  are  computed  at  maximum 
depth  and  near  the  fin  surface. 

iv.  While  cracks  emanating  from  the  fin  tip  on  the 
axis  of  symmetry  remain  essentially  planar,  the 
off-axis  cracks  generate  curved  surfaces. 

The  preliminary  study  described  above  suggests 
that  the  use  of  off-axis  cracks  may  be  an  unnecessary 
complication  in  predicting  crack  growth  and  SIF  de¬ 
termination  due  to  the  substantial  irregularities  and 


highly  elevated  loads  required  to  induce  sufficient  crack 
growth  in  off-axis  cracks  to  shed  the  shear  mode  for 
the  geometry  studied.  However,  further  experiments 
are  underway  in  order  to  further  extend  the  results. 
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APPENDIX  A-  Frozen  Stress  Photoelasticity 

When  a  transparent  model  is  placed  in  a  circularly 
polarized  monochromatic  light  field  and  loaded,  dark 
fringes  will  appear  which  are  proportional  to  the  ap¬ 
plied  load.  These  fringes  are  called  stress  fringes  or 
isochromatics,  and  the  magnitude  of  the  maximum  in¬ 
plane  shear  stress  is  a  constant  along  a  given  fringe. 

Some  transparent  materials  exhibit  mechanical 
diphase  characteristics  above  a  certain  temperature, 
called  the  critical  temperature  ( Tc ).  The  material, 
while  still  perfectly  elasticjwill  exhibit  a  fringe  sensitiv¬ 
ity  of  about  twenty  times  the  value  obtained  at  room 
temperature,  and  its  modulus  of  elasticity  will  be  re¬ 


used  throughout  this  paper. 

APPENDIX  B 

MODE  I  ALGORITHM 

Beginning  with  the  Griffith-Irwin  Equations,  we  may 
write,  for  Mode  I,  for  the  homogeneous  case, 

=  7~-rfij(Q)  +  v-j  ( i,j  =  n,z )  (1) 

{/nr )  2 

where: 

&ij  are  components  of  stress 
is  SIF 

r,0  are  measured  from  crack  tip  (Fig.  B-l) 


duced  to  about  one  six-hundredth  of  its  room  temper¬ 
ature  value.  By  raising  the  model  temperature  above 
Tc,  loading,  and  then  cooling  slowly  to  room  temper¬ 
ature,  the  stress  fringes  associated  with  Tc  will  be  re¬ 
tained  when  the  material  is  returned  to  room  temper¬ 
ature.  Since  the  material  is  so  much  more  sensitive  to 
fringe  generation  above  Tc  than  at  room  temperature, 
fringe  recovery  at  room  temperature  upon  unloading 
is  negligible.  The  model  may  then  be  sliced  without 
disturbing  the  “frozen  in”  fringe  pattern  and  analyzed 
as  a  two-dimensional  model  but  containing  the  three- 
dimensional  effects.  In  the  use  of  the  method  to  make 
measurements  near  crack  tips,  due  to  the  need  to  re¬ 
duce  loads  above  critical  temperature  to  preclude  large 
local  deformations,  and  the  use  of  thin  slices,  few  stress 
fringes  are  available  by  standard  procedures.  To  over¬ 
come  this  obstacle,  a  refined  polariscope  is  employed  to 
allow  the  tandem  use  of  the  Post  and  Tardy  methods 
to  increase  the  number  of  fringes  available  locally. 


<7?.  are  non-singular  stress  components. 

Then,  along  9  =  7t/2,  after  truncating  aij 

rmax  _  _JSl_  .  T°  _  KAP  (2) 

nz  (8tt  r)§  (87rr)i  U 

where:  r°  ==  /(eq?)  and  is  constant  over  the  data  range 
Kap  =  apparent  SIF 
r™ax  —  maximum  shear  stress  in  nz  plane 

•  kap  _  K1  ,  VSt0  /n  i  /,•) 

a(i ra)i  ro)  2  9  \a) 


where  (Fig.  B-l)  a  =  crack  length,  and  a  =  remote 
normal  stress 


Kap  fr  . 

1,e*  T7 — nr  vs*  \  ~  1S  lmear. 
<7(7Tg)2  V  a 


Since  from  the  Stress- Optic  Law 

7_max  —  — ^  where 
‘nz  21  W1ACAC 

n  =  stress  fringe  order 


In  fringe  photographs,  integral  fringes  are  dark  in  a 
dark  field  and  bright  in  a  bright  field.  Bright  fields  are 


/  =  material  fringe  value 
t  =  specimen  (or  slice)  thickness 
and  from  Eq.  2 

Kap  =  r“ax(87rr)*  =  ^(8tt  r)i, 

then  Kap  (through  a  measure  of  n)  and  r  becomes 
the  measured  quantity  from  the  stress  fringe  pattern 
at  different  points  in  the  pattern. 

In  the  present  study,  instead  of  normalizing  K 
with  respect  to  a(na)1/2,  we  have  selected  p^/na/Q  as 
the  normalizing  factor  where  y/Q  is  an  elliptic  integral 
of  the  second  kind  approximated  here^as  shown  in 
Table  I.  An  example  of  the  determination  of  i*i  in 
Table  I  from  test  data  is  given  in  Fig.  B-2. 

MIXED  MODE  ALGORITHM 

The  mixed  mode  algorithm  was  developed  (see 
Fig.  B-3)  by  requiring  that: 

Urn  { {Snrm)1/2 (Ku A2 , rm , 0m , ) |  =0 
©mm-^©L  ^ 

(4) 

which  leads  to: 

(5) 

By  measuring  0^  which  is  approximately  in  the  direc¬ 
tion  of  the  applied  load,  K2/K1  can  be  determined. 


Knowing  K*  and  0^,  K\  8z  K<i  can  be  deter¬ 
mined  from  Eqs.  5  and  6.  Details  are  found  in  Smith 
and  Kobayashi  (1993). 


Fig.  B-l  Near  Tip  Notation  for  Mode  I 


Model  3,  short  symmetric  crack, 
a/c  =  0.5,  a  ft  =  0.41,  depth  (d) 


Fig.  B-2  Determination  of  F\  from  Test  Data 


Then  writing  the  stress  optic  law  as: 

—max  -  LH  -  ^AP 
2t  (&irr)i 


one  may  plot 


^AP^P-  vs.  Jr  fa  as  before,  locate  a 
pjna 

linear  zone  and  extrapolate  to  r  =  0  to  obtain  K*. 
Now  K *,  may  be  written  as: 


K*  =  [(K^ine^+lKicosQlf+iKisinQl)2}*  (6) 


ge  order  (N) 


Fig.  4  Fringe  Patterns  Near  Critical  Loci  at  Fin  Tip 


fringe  order  estimation  along  the  axis  of  symmetry 
and  offraxis  of  symmetry  around  the  fin  tip 
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Fig.  5  Convergence  of  Fringe  Orders  at  Critical  Loci  in  Uncracked  Fin 
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CFS -corner  of  fin  surface 
C- crack  front 


I  -  initial  crack 
F  -  final  crack 


FS-fin  surface 


MF  =  4.46 


Fig.  7  A  Long  Off-Axis  Crack  (Model  2) 


